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I. INTRODUCTION 



Entanglement, the characteristic trait of quantum mechanics according to Schrodinger, has been identified as a 
resource central to many of quantum information processing tasks There have been tremendous progress along 
the characterization and the quantification of entanglement using various methods in both bipatite and multipartite 
settings [3| . Entanglement has also been studied in many-body systems and connection to quantum phase transitions 
has also been explored [3j . The notion of entanglement has also found its application in density matrix of renormal- 
ization groups (DMRG) and has enabled much recent progress in numerical techniques of DMRG in low dimensions 
and dynamics |4!|. 

Despite these advancements, the characterization and the quantification of entanglement are far from complete. For 
example, even for multipartite pure states, it has been a long standing question whether there exists a finite minimal 
reversible entanglement generating set (MREGS) @. The existence of MREGS would enable the generalization 
of the entanglement of distillation Jj] and formation 0, § to multipartite systems and would also provide a 
better characterization of multipartite entanglement. Given that the standard measures of entanglement, such as the 
entanglement of formation 0, [1] and distillation 0, Q have not be properly generalized to multipartite settings, the 
study of quantifying multipartite entanglement via other measures is indispensable. 

The relative entropy of entanglement (REE), introduced by Vedral et al. |Tol.lTl|. provides an alternative measure of 
entanglement. In bipartite settings, it is shown to be a lower bound for the entanglement of formation and an upper 
bound for the entanglement of distillation. Its regularized version has recently been shown to possess connection to 
the second law of thermodynamics fT5]. Moreover, REE applies straightforwardl y to multipartite settings, and thus 
its study in the multipartite settings is important. The calculation of REE (lol. Illl| by analytic means, however, is 
still a challenging task, even for pure states, and is often relied on numerical computation. Despite the progress made 
in two-qubit systems by Ishizaka [13], an analytic formula for REE is still elusive. 

In this paper, we shall focus on the study of REE, complemented by two other measures. These other measures, the 
geometric measure of entanglement (based on the geometry of Hilbert space [H, U^, and also based on the Grover 
search [13, [Si i hence, also known as the Groverian entanglement) and the generalized robustness of entanglement [l^ , 
besides providing different perspectives of entanglement on their own, can also be used to provide connections to REE, 
or more precisely, lower and upper bounds, respectively. The main results of this paper are the establishment of the 
proof for conjectures for analytic expressions of REE for several families of multipartite mixed states, thus providing 
nontrivial examples where REE is obtained analytically. 

The structure of the present paper is as follows. In Sec. |TT] we review the three entanglement measures considered 
in the paper: the relative entropy of entanglement, the geometric measure of entanglement, and the robutsness of 
entanglement. We explore connections among the three, in both pure- and mixed-state settings in Sec. lIIII and Sec. lIVi 
respectively. In Sec.|V]we prove the conjecture for multi-qubits (20| and we also generalize the result to the multi-qudit 
setting. With these analytic results, we discuss some applications in Sec. IVIII In Sec. IVIIII we prove the analytic 
formula of REE for another family of multi-qubit mixed states [2l|, which include Diir's bound entangled states [23 |. 
In Sec. llXI we give some concluding remarks. 
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II. ENTANGLEMENT MEASURES 

In this section we briefly review the three measures considered in the present paper: the relative entropy of entangle- 
ment, the general (global) robustness of entanglement and the geometric measure of entanglement. 

A. Relative entropy of entanglement 

The relative entropy ^(/cilItT) between two states p and tr is defined via 

5(p||a) = Tr(plog2p-plog2a), (1) 

which is evidently not symmetric under exchange of p and a, and is non-negative, i.e., S'(p||(t) > 0. Note that the log 
function is base-2 throughout this paper. The relative entropy of entanglement (RE) for a mixed state p is defined to 
be the minimal relative entropy of p over the set of separable mixed states [13, [Hi : 

Ek{p) = minS'(p||CT) = min Tr (p logj p - p log2 <t) , (2) 

where V denotes the set of all separable states. 

The regularized relative entropy of entanglement is defined as 

E^ip) EE hm -i?R(p®"). (3) 

It is shown that in the bipartite settings the regularized relative entropy of entanglement plays an analogous role to 
the entropy in thermodynamics 12]. However, the calculation of the regularized relative entropy of entanglement is, 
in general, much more difficult than for the non-regularized case. 

In general, the task of finding the REE for arbitrary states p involves a minimization over all separable states, and 
this renders the computation of the REE very difficult. For bipartite pure states, the REE is equal to the entanglement 
of formation and of distillation. But, despite recent progress [l^l, for mixed states — even in the simplest setting of 
two qubits — no analog of Wootters' formula [1] for the entanglement of formation has been found. Things are even 
worse in multipartite settings. Even for pure states, there has not been a systematic method for calculating their 
relative entropy of entanglement. It is thus worthwhile seeking cases in which one can explicitly obtain an expression 
for the REE. 

An alternative definition of RE is to replace the set of separable states V by the set of postive partial transpose 
(PPT) states Vpp*: 

VPP' = {a\\a^ =a, a> 0, Tr(CT) = 1, a^^ > 0}, (4) 

where PT denotes partial transpose with respect to any bi-partition of parties. The REE thus defined, as well as 
its regularized version, gives a tighter bound on distillable entanglement. There has been important progress in 
calculating the RE (and its regularized version) with respect to PPT states for certain bipartite mixed states; see 
Refs. [235 for more detailed discussions. For multipartite settings one could also use this definition, and define the set 
of states to optimize over to be the set of states that are PPT with respect to certain or all bipartite partitionings. 
However, we shall use the first definition, i.e., optimization over the set of completely separable states, throughout 
the discussion of the present paper. 

B. General robustness of entanglement 

The general robustness of entanglement [l^ is a measure of how sensitive the entanglement is to mixture of states. 
It is defined as 

R{p)EEmint (5) 
such that there exists a state A so as to render the following state separable: 

1 



:{p + tA). (6) 



1 + t' 

The logarithmic robustness of p is defined as 

LR{p) = log, (l + R{p)). (7) 
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C. Geometric measure of entanglement 



We continue by briefly reviewing the formulation of the geometric measure in both pure-state and mixed-state settings. 
Let us start with a multipartite system comprising n parts, each of which can have a distinct Hilbert space. Consider 
a general n-partite pure state (expanded in the local bases {Icp.'}): 



E 

We can compare this state to the set of general separable pure state 



1^-)= E Xp.P.-pM\^eg^---e(-J). (8) 



|0)^^Sj^«)=^S^(^c«|e«)), (9) 
Pi 

and define the maximal overlap of \ip) with the closest product states as follows, 

A„,ax(^) =niax|(0|^)|, (10) 
•P 

where |0) is an arbitrary separable pure state defined above. In Ref. [lB|, the particular form £^[^2 = 1 — A^j^^(V') = 
sin^ ^niin was defined to be the geometric measure of entanglement (GME) for any pure state Here, we shall be 
concerned with the related quantity Eiog{ip) = — 2 log2 Aniax(^), which we shall show to be related to the two other 
measures. 

We remark that an alternative perspective from the Grover search also leads to the geometric measure [TtI [isj , and 
it is also known as the Groverian entanglement. Furthermore, a hierarchy of entanglement can be obtained if one 
consider the separable states to be product states among an appropriate partition of particles into fc-parties [l^, [lB| ; 
see recent works by Shimoni and Biham and by Blasone et al. [1^. But throughout this paper we shall focus on 
the completely product states, and the separable mixed states will refer to n-separable states. 

Given the definition of entanglement for pure states just formulated, the extension to mixed states p can be built 
upon pure states via the convex hull construction (indicated by "co"), as was done for the entanglement of formation; 
see Ref. Q. The essence is a minimization over all decompositions p = J2iPi into pure states: 

E{p) = {coEpa,.c){p) = min V impure (IV'i))- (11) 
{pi.V-i} — ^ 

This convex hull construction ensures that the measure gives zero for separable states; however, in general it also 
complicates the task of determining mixed-state entanglement. The specific form that we are concerned in this paper 
is the following 

Eiogip)= min V [ - 2 log Amax('0i)] ■ (12) 
{pi,^i} 

We remark that an alternative way to define the mixed-state entanglement measure via purification have been 
developed by Shapira et al. [26| , continuing along the idea of maximizing the Groverian search outcome. But in the 
present paper, we shall focus mainly on the definition by Eq. (|12p . 

Illustrative examples: We examine some pure-state examples, whose mixture will be considered later in the paper. 
First, one can classify permutation-invariant pure states, as follows: 

|5(n,fc)) EE S\0_j_01_--1) = ^Y.^^\0_■3L:J)^ (13) 



^ k n-k * k n-k 

where S represents symmetrization under permutations, 11.; denotes arbitrary permutation of n objects, and CJ} = 
n\/k\(n — fc)!. Intuitively, as the amplitudes are all positive, one can assume that the closest separable (equivalently, 
Hartree) state is of the form (which is rigorously proved in Ref. [27l |) 

10) ^(VP|0) + 71^11))^", (14) 
for which the maximal overlap (w.r.t. p) gives the entanglement eigenvalue for |S(n, fc)): 
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More generally, for n parties each of which being a d-level system, the state 



\S{n-k)) = V ''^'^,''^''' 'g|l_^2_^...d_^) (16) 

ki k-2 fcd 



has the entanglement eigenvalue 



Although the above states were discussed in terms of the GME [l6[, we shall, in the following section, show the 
rather surprising fact that the LR and REE of these example states, are given by the corresponding expression: 

-21og2 Amax- 

One of the main results of the present paper concerns the two following families of mixture of symmetric states: 



pm) = (18) 

fc=0 

= Y.n\S{n,k)){S{n,k)\- (19) 



III. CONNECTION BETWEEN THE THREE MEASURES: PURE STATES 

In bipartite systems, due to the existence of Schmidt decompositions, the relative entropy of entanglement of a 
pure state is simply the von Neumann entropy of its reduced density matrix. However, for multipartite systems there 
is, in general, no such decomposition, and how to calculate the relative entropy of entanglement for an arbitrary pure 
state remains an open question. We now connect the three measures by inequalities. 

A. Relative entropy of entanglement and geometric measure of entanglement 

Let us begin with the following inequality: 

Inequality 1 [20l. [28t: For any pure state It/j) with entanglement eigenvalue Aniax(V') the quantity — 2 log2 Amax(V') is 
a lower bound on the relative entropy of entanglement of j?/'), i.e., 

£^r(|^)(V'I) > -21og2A^ax(^). (20) 

Proof: From the definition ^ of the relative entropy of entanglement we have, for a pure state j^/"), 

£^r(|^)(^|) =min-(V'|log2a|V') =-max(^|log2a|V'). (21) 

Using the concavity of the log function, we have 

(V|log2a|V) <log2((Vkl^)) (22) 

and, furthermore, 



We then conclude that 



max('0| log2 cr|?A) < maxlog2(('(/'k|'0))- (23) 

E^mm > -maxlog2((Vk|^)). (24) 
As any a €V can be expanded as a = ^iPi\(t>i){(l>i\, where |(/>?)'s are separable pure states, one has 

= ^K|('/'d^>l' < ALxW, (25) 

i 

and hence we arrive at the sought result 

£^r(|^)(V'I) > -21og2A^ax(^). (26) 

□ 
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B. Geometric measure and generalized robustness 

Cavalcanti has obtained the resuh that robustness of entanglement can be shown to provide an upper for the 
geometric measure j29| . The insight he provided is the connection of both measures to entanglement witness. We 
shall repeat his result for pure states here. The connection of the generalized robustness to entanglement witness is 
the result of Brandao [s^l : 

R{p) = max{0, - minTr(VFp)}, (27) 
w 

where M < 1 and T'cipsEp) > 0. Here we shall be concerned with pure states, and hence 

i?(|^A)(^|)=max{0,-min(V'|M^|^))}. (28) 

The connection of the geometric measure of entanglement to entanglement witness is as follows. The operator 

W={XH~-\yj){^p\)/X^ (29) 

satisfies the condition Tr{WpsEp) > as long as > A'^^^{^|J) Using this and take W* = (Aniax(^/')^1 — 

|V')(V'l)/Amax(V')i it will provide a lower bound on R{\%p) {tp\) , and thus 

Rmm > ^77^%?- (30) 



Aiax(V') 



Equivalently, we have the following 
Inequality 2 [2^1 : 



M(|V)(^|)>-log2ALx(V'). (31) 
It turns out that a stronger lower bound on the robustness can be obtained using the relative entropy of entanglement, 



first proved by Hayashi et al. 3l|, which we now describe 



C. Generalized robustness and relative entropy of entanglement 

We repeat here the proof that for any pure state (with p^, = |?/')(V'|), 
Inequality 3 fsij: 

Eni\i;){i;\)<LR{\^P){4'\). (32) 
Proof. Suppose R{p^) — tm- There exists a state A such that the following state is separable 

<JS = -^(\^){^\+Un^)- (33) 

Thus we have 



ERm{^)<s{\^){^\ 



o-s = -(V'llog- 



< -(V.|logMM|^)=-log-^ = Li?(|^)(V|). (34) 



□ 



Inequalities 1 to 3 summarize to 
Inequality 3': 



Li?(|V)(^|) > EM{4'\) > -21og2 A^ax(V')- (35) 

The degree of difficulty of calculating the corresponding measures becomes higher from the r.h.s. to the l.h.s. Moreover, 
upper bounds on LR and Sr, are easy to obtain, as they are defined via minimization. If one can construct an upper 
bound on these measures such that it matches — 21og2 Amax(V')j the inequalities above become equalities. This is how 
the analytic expression of left two measures have been obtained for symmetric and antisymmetric states as well as 
some graph states [13, [13, [lH, [12] . A very thorough analysis using group theory on whether the inequalities become 
inequality has been provided by Hayashi et al. [27|. 
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D. Illustrative examples 

We now examine illustrative states in the light of the above discussions, thus obtaining the expression of the three 
different entanglement measures. We begin with the permutation-invariant states \S(n,k)) of Eq. (|13p . for which 
Amax was given in Eq. (jlSp. The above theorem guarantees that E-[i{\S{n, fc))) > — 21og2 A„iax("': k). To find an upper 
bound we construct a separable mixed state 

^m))m)l (36a) 
1^(0)) = (vp|0)+e'*yT^|l))®", (36b) 



with p chosen to maximize k))\ \ = \/C^p^{^—py^ which gives p = k/n. Direct evaluation then gives 

n 

=E^JV(l-p)'"-^'l^(n,j))(5(n,j)| = ALxKfc)|^Kfc))(^(n,fc)|+r^, (37) 
i=o 

and S'(/9||ct) = — 2 log2 A,„ax('T-, fc), where p — |S'(n, fc)) (5(n, fc)|, A,„ax('T-,fc) is given in Eq. (|15p . and we have used 
to denote the those terms with (j ^ k). The upper and lower bounds on E'r coincide, and hence we have that 

£;R(|5(n,fc))) =-21og2A,„ax(n,fc). (38) 

Moreover, from the perspective of the generalized robustness, we immediately have an upper bound i,„ on the 
robustness, 

^ =ALxKfc)- (39) 



1 + 1... 

Hence, we obtain an upper bound on LR, 

LR{n, k) < log2(l + t,„) = -2 log2 A,nax(n, k). (40) 

This then gives equality to all three measures for the symmetric state \S{n,k)). The same consideration also holds 
for \S{n,k)), namely, 

LR{n,k) = Enin,k) = -2 log2 A„iax(n, ^). (41) 

IV. CONNECTION AMONG THE THREE MEASURES: MIXED STATES 

Cavalcanti and Hayashi et al. have provided inequalities relating these three measures for general mixed states. In 
their inequalities, the geometric measure is generalized to 

G(p) = -log2{maxTr(H}. (42) 

We shall see below that inequalities in terms of -Eiog can also be derived and shall provide an alternative inequality. 

Let us first review the inequality shown by Cavalcanti [2^ : 
Inequality 4 [2^: 

LR(p) > G{p). (43) 
Recall that we take V to be the set of separable states of the form 

Y.P^\^^)^^^\®\^f)^^?\®---^ (44) 



and hence G{p) can be simplified to be 



G(p) = - log2{ max (0|p|<^)}, (45) 

(pGProd 
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where Prod denotes completely product states. 

Proof: Again, we use the relation ([77]) for the robustness and the entanglement witness obtained by Brandao [30| . 
For any state p, construct a witness 

Wx = t- ^p, (46) 

where > A^^^^ = maxo-eP Tr (per), in order for the condition Tr{W\a) > to be satisfied. Plug in the above 
expression for witness into Eq. (|27p . we obtain 



1 

A^ 



R{p) > -MWx^J = -l + -^. (47) 



max 



This is equivalent to 

Li?(p) >-log2ALx = G(p). (48) 



□ 



The inequality derived by Hayashi et al. is as follows, 
Inequality 5 [31i] : 



log^r{p)>Enip) + S{p)>G{p), (49) 
where, instead of the robustness of p itself, the robustness of the support of p is considered: 

r{p)^\Pp\[l + R{P,/\Pp\)l (50) 

where Pp is the support of p, namely 

Pp- ^ |A,)(A,|, (51) 

A, 7^0 

where p has the spectral decomposition p — Xi\Xi){Xi\, and \Pp\ = Tr(Pp). 

We now slightly extend the first part of the inequality. By the definition of robustness, there exist a state A such 
that the state 

Lo = (52) 

l + t ^ ' 

is a separable state, where t = R{p) for convenience. This gives 

Er{p) + S{p) = min-Tr(plog(T) 

< -Tr(plogtj) 

p + tA 
l + t 
P 



= -Tr yplog 
< -Tr (plo, 



1 + */ 

= S{p) + LR{p). (53) 

This gives 

LR{p) > Er{p), (54) 

where the l.h.s. is given by the logarithmic robustness of p, rather than that of its support Pp/\Pp\. The second 
inequality of Hayashi et al. is elementary to prove: 

Er{p) + S{p) ~ — max Tr(ploga) > -log{ maxTr(p(T)}, (55) 
where we have used Tr(p log cr) < logTr(ptT). 
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We can also provide an alternative of the second-part of Inequality 5: 

Er{p) > E,oM - S{p). (56) 

Proof: 

Suppose p — is the optimal decomposition for £'iog(p), namely, 

^^iog(p) = -Y.P- l°g2 ALx(V'»)- (57) 

i 

Using the definition of REE, we have 

Er{p) + S{p) = -maxTr(plog(T) > - max Vp,Tr(|V'»)(^»| loga.). (58) 

i 

Using 



we have 



Summing up we have proved 
Inequality 6: 



maxTr(|V'^)(V'^|logc^,) = logA^,^,,(V'^), (59) 



Er{p) + S{p) > -^paogAL.(^») =i?iog(p). (60) 

i 

□ 

LR{p) > Er{p) > - S{p). (61) 



Let us compare Inequality 6 to Inequality 5, especially the first parts. Viewed as an upper bound on E^, the 
Inequality 5 can be rewritten as 

\og^r{p)-S{p)>En{p). (62) 
It can happen that the r.h.s. is zero whereas the l.h.s. is still nonzero, as exemplified for the two-qubit Werner state, 

Pw(7)=7lV'")(V'"l + ^^l4x4, (63) 

where \ip~) is the two-qubit singlet state. The state is entangled for 7 > 1/3, and unentangled for 7 < 1/3. For 
< 7 < 1, the support of p\Y is Pw = l4x4: and the corresponding state is a completely mixed state, hence, 
possessing no entanglement. This leads to r{pw) — 4, and 

^ ^^ 1 + 37, 1 + 37 3(1-7), 1-7 

At 7 = 1/3, Enipw) = 0, but log2r(p) - S{p) = 1 - log2(3)/2 k. 0.208. Although LR{p) is generally not easy to 
calculate, it becomes zero when p becomes separable. For the second-part of Inequality 6, we remark that for the 
states in Eqs. (fT8|) . (fT9|) and ()144p the lower bound can be tightened to 

Er{p) > E.o^ip), (65) 

which will be shown later. Therefore, Inequality 6 slightly extends previous results by Cavalcanti [29*1 and Hayashi et 

ai. mi. 

V. RELATIVE ENTROPY OF ENTANGLEMENT FOR MIXTURE OF SYMMETRIC STATES 

A. Multi-qubits 

In Ref. [16] the procedure was given to find the geometric measure of entanglement, i?sin^; ^or the mixed state 
comprising symmetric states: 

p(M) = E^"^-l^(^'^))(^("'^)l- (66) 
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Here, we focus instead on the quantity Eiog, but the basic procedure is the same. The key point is to find the 
entanglement eigenvalue A„({g}) for the pure state 

5^v^|5(n,fc)), (67) 

k 

thus arriving at the quantity 

f({g}) = -21og2A„({g}). (68) 
Then the quantity E\c,g for the mixed state ([66]) is actually the convex hull of the expression (|68|) : 

£^iog,(p(M))-co£:(M). (69) 

In Ref. ^2(1, an attempt to calculate the REE for the states (|66p was made and a conjecture for the REE was made. 
In this section, we review the construction of the conjecture and prove it to be correct. Let us consider the state 
formed by mixing the separable pure states 1^(6*, </>)): 

<j{e) = / g|C(0,0))(e(^,</>)| =X^Q'cos2^0sin2("-'^-)0|5(n,fc))(5(n,fc)|, (70) 

k=Q 

where 

|e(0,0)) = (cos0|O)+e^"*sin0|l))®". (71) 
Allowing 6 to vary, we then minimize the relative entropy between p{{p}) and (j{9), 

s ipi{p})Me)) = J2pk log l\ (72) 

^ cos2'= sm"^^" 



with respect to 0. We arrive at the stationarity condition 



Due to the convexity of the relative entropy, namely, 



\ i i / i 



(73) 



(74) 



we can further tighten the expression of the relative entropy by taking its convex hull. The convexification process 
also results in the corresponding separable state, i.e., the knowledge of the coefficients r^'s in a* 

a* =^r,|5(n,fc))(^(n,A:)|. (75) 
fe 

The arrived upper bound on REE for the mixed state p{{p}) was conjectured (shall be proved below) in Ref. [l^l to 
be the exact REE: 
Theorem 1: 

En {p{{p}))^ CO F{{p}), (76) 

where 

■F'(M) = y'Pfel0g2 73 9t 2(n-k) /, ^ l^Pfc ^Qg2 r^n k^ " \n-k ' ('"^) 

^ C^'cos^'^^esin^^" ^ CJ!a'^(n-ap *= 

where the angle satisfies Eq. ([75|). CJ^ =n\/ [kl{n — A:)!), and a = J2kPk k- 
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B. Proof of the theorem 

1. Symmetry considerations 

We shall discuss the symmetries possessed by the states p({p}), and these shall reduce the set of separable states 
that we need to consider. We begin by noting that the states p{{p\) are invariant under the projection 

^^■■P^ J ^U{cf,ryU{cf>)^^- (78) 

with C/(0){|O), |1>} ^ {|0),e-''^|l)}. VoUbrecht and Werner [H] have shown that in order to find the closest separable 
mixed state for a state that is invariant under projections such as P, it is only necessary to search within the separable 
states that are also invariant under the projection. 

We can further reduce the set of separable states to be searched by invoking another symmetry property possessed 
by p{{p}): these states are also, by construction, invariant under permutations of all parties. Let us denote by 11^ one 
of the permutations of parties, and by Iii{p) the state obtained from p by permuting the parties under 11^. We now 
show that the set of separable states to be searched can be reduced to the separable states that are invariant under 
the permutations. To see this, suppose that p is a mixed state in the family (|66p . and that a* is one of the closest 
separable states to p, i.e.. 



As p is invariant under all 11^, we have 



EK{p)^rmnS{p\\a)^Sip\\a*). (79) 



^R('°) = ]^E^(H|n.(^*)), (80) 



where iVn — nl is the number of permutations. By using the convexity of the relative entropy we have 

Enip) > S (^p\\[J2^^i<J*)/Nu]^ . (81) 

However, because of the extremal property, Eq. ([79)l . the inequality must be saturated, as the left-hand side is already 
minimal. This shows that the state under the projection 

P2:^*-]^En4^*) (82) 

i 

is also a closest separable mixed state to p, and is manifestly invariant under all permutations (noting that P2 preserves 
separability). Thus, we only need to search within this restricted family of separable states, namely separable states 
invariant under Pi and P2. 

In fact, one can use the group representation theory (see, e.g. Ref. [13, [13]) to characterize all states that invariant 
under Pi and P2. It turns out that any such state (invariant under both Pi and P2) can be written as mixture of 
all \S{n,k)) and those other basis states belonging to other irreducible representations for the group Sn- However, 
we have not been able to prove our theorem using only symmetry argument, for it is not trivial to characterize 
all separable states of this form. But we can prove it by an algebraic approach already developed by Vedral and 
Plenio [li|. 

It is not difficult to see that the set Vs of all separable mixed states that are diagonal in the basis of {|S'(n, fc))} 
(basis states for the totally symmetric subspace Vs) can be constructed from a convex mixture of separable states in 
Eq. ([70]) . That is, for any (Tg S I?s we have a decomposition 

a, =^i,c7(&0, (83) 

i 

where ti > 0, J2i ti — ^i a-nd <j{9i) is of the form ((70| . This is because the separability of the states (|66p implies that 
there exists a decomposition into pure states such that each pure state is a separable state. Furthermore, because 
{[^(n, fc))} are eigenstates of p{{p}), the most general form of the pure state in its decomposition is 



^y^e'*'=|5(n,fc)). (84) 
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This pure state is separable if and only if it is of the form (j7ip . up to an overall irrelevant phase. As p{{p}) is invariant 
under the projection Pi ([75)) . a pure state in Eq. (|7ip will be projected to the mixed state in Eq. ([70)1 under Pi. Thus, 
every separable state that is diagonal in {|S'(n, k))} basis can be expressed in the form ((83|) . 

Hence, our construction of a* (via any necessary convexification) ensures that it achieves at least the minimum (of 
the relative entropy) when the separable mixed states are restricted to Vs. However, in order to prove the conjecture, 
one would still need to show that the expression is also the minimum when the restirction to Vs is relaxed to the set 
of separable states invariant under Pi and P2, unless we can employ further argument, symmetry or not, to reduce 
to T>s the set of separable states that we need to consider. 



2. The proof 

The symmetry considerations presented in the previous section only reduce the set of separable states to a smaller 
one possessing symmetry. In particular, any state that is invariant under (1) any permutation of the n parties and 
(2) projection under average of local phase as in Eq. l|78|) can be written as a mixture of the basis states in all the 
irreducible representations of U(2). However, the symmetry alone has not led us to the proof of our theorem. Here, 
we provide a more direct algebraic proof. 

Since we have shown that a* gives minimum of S'(/3||(t) when the separable states are restricted to being diagonal in 
\S{n, fc))'s, we now show that it is indeed a local minimum when we lift that restriction. According to the discussions 
by Vedral and Plenio [Tlj . this means that we need to show that adding any separable state as'. cr{x) = {l — x)a* +xas, 
the quantity /, 

fix,as) = Sip\\{l-x)a* +xcTs), (85) 
has a local minimum at x = 0. This is equivalent to show that 

|^(0,a.)>0. (86) 

It is straightforward to calculate the l.h.s. of the above expression, as was done in Ref. [ll|, which gives 

1^(0, a,) = 1 - dtTT[pia* + i)- V,(a* + t)-'] . (87) 

Taking a* as in Eq. ([75]) and p in Eq. we have 

1^(0, a,) = 1 - P^{S{n, k)\as\S(n, k)). (88) 

As a* is constructed to be the minimum in the restricted (i.e., the totally symmetric) subspace, this means that for 
|$s) — (cos 6*10) + sin6'|l))'*" in this subspace, we automatically have by construction 

^(0,|$,)($,|) = l-V — C^cos2^0sin2("-^)0> 0. (89) 

ox ^ Tfe 

K 



As 9/(0, a)/dx is linear in ct, namely. 



g(0,5:fta.) = E^'»^(0'^^)' (90) 



in order to show that df{Q,cr)/dx > holds for arbitrary separable state cr, it is sufficient to show that it holds for 
arbitrary pure separable state a = |$)($|, which is what we are about to do. To be more explicit, we shall show that 
for any arbitrary separable pure state 

l'i>> = ®■=l(^/^|0) + VT^e^^^ll)), (91) 

we also have 

|^(0,|$)($|) = l-^^(5(n,fc)|$)($|5(n,fc)) >0. (92) 
""^ k=0 
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0.6 0.8 



FIG. 1: Comparison of F (solid curve), E-r. = co F (partly dashed line, when convexification is necessary) and the numerical 
value of En (dots) for the states p3;o,i(s), p3;o,2(s), and P3;i,2(s) (from top to bottom). 



Note that |$) is the most general separable multi-qubit pure state (up to an irrelevant global phase). 
The key point is then to evaluate 



$fc = (5(n,fc)|$). 



Note that 



This gives 



HieSn 



k n—k 



(93) 

(94) 



$fc = (5(n, fc)|$) = ^ E • • • \/to;wVl-9n,(fc+i)e''^°^"=+^' . . . ^1 - 9n,(„)e''^"*w. (95) 



Thus, 



and 



Hi 



< 7-^ [E«n.(i) • • •ta.Cfe)] [E(l - 9n,(fe+i)) • • • (1 - 5n, 

^ '-' IT. TT. 



(96) 



(97) 
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where we have used the Cauchy-Schwarz inequahty. Using the Maclaurin inequahty (38|: 

fe 



-^xn.a)...xn.(.)<(^^^j ' 



Hi 

we have that 

l^'fcl^ < Ck{qf{l - q)'^-^ = CI cos^'^ 6isin2("-'=) 9, (99) 

for cos^ 9 = q. This means that 

^(0,|$)($|) = l-V ^|$fep > l-V ^Q'cos2'=0sin2("-'^-)e>O, (100) 

and that cr* is indeed the closest separable state to p in Eq. (|66p . Hence, Theorem 1 is proved. 

C. Examples 

We illustrate the established expression of REE for the state p({p}), making the restriction to mixtures of two 
distinct n-qubit states \S{n,ki)) and \S{n,k2)) (with fci ^ k2)'. 

p^,k,.M (s) = s\S{n, ki)){S{n, fci)| + (1 - s)\S{n, k2)){S{n, k2)\. (101) 

We first investigate the two-qubit (i.e. n — 2) case. Besides the trivial mixture, P2;0,2, there is only one i neq uivalent 
mixture, P2;0,i{s) [which is equivalent to P2;2,i(s)], which is the so-called maximally entangled mixed state [36L[37j (for 
a certain range of s) 

P2;0,l = S |11)(11| + (1 - .S)|V|/ + )(M/+|, (102) 

where 1*+) = (|01) + 10)/V2 is one of the four Bell states. The function F for this state [denoted by i^2;0,i(s)] is 

F2;oAs) = S \0g2 + {1-S) l0g2 ^ , (103) 

which is convex in s. Hence, it is exactly the expression for the relative entropy of entanglement for the state P2:0,i 
found by Vedral and Plenio [lH|. 

For n = 3 there are three nontrivial inequivalent mixtures: P3;o,i(s) [equivalent to P3-3.2{s)], P3;o,2(s) [to /53;3,i(s)], 
and P3;i,2(s) [to P3;2,i(s)]. In Fig.[l]we compare the function F in Eq. ([77)) . its convex hull coF, and numerical values 
of £'r obtained using the general scheme described in Ref. [ll[ extended beyond the two-qubit case. 

For n = 4 there are five inequivalent nontrivial mixtures: P4;0,i(s), P4;0.2(s), Pi-.o.iis), /34;i.2(s), and p4-i.3{s). In 
Figs. [3] and m we again compare the function F in Eq. ([77)1 . its convex hull coi^ = E^i, and numerical values of E'r. 

For the states that we have just considered, we now explicitly give the formulas for E^i suggested by the theorem. 
For the three-qubit mixed state P3;2,i('S)7 its -Br is 



9s 9(1 — s) 

(l + s)2(2-s)+^^"')^°S2 (2 _,)2 (! + ,)• 



For P3:0,i(s), it is 



For P4;0,i it is 



For p4;i,2(s), it is 



27s 9 

^^°s^(2T^ + (^-^)^°s^(2T^- 



256s 64 , 

^^°^^(3T^ + ('-^^^°^^(3T^- 



64s , , 128(1 - s) , , 

(2-s)(2 + s)3 + (1-^) 3(2-s)^(2 + s)^ - ^^^^^^ 
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FIG. 2: Comparison of F (solid), its convex hull, i.e., E^, and the numerical value of _Er for the states p4;o,i(s), P4;i,2(s), and 
P4;i,3(s) (from top to bottom). In these cases, E-r = F. 



For p4;i,3(s), it is 



Sl0g2 



64s 



(3-2s)(l + 2s)3 



+ (l-s) log2 



64(1 - s) 
(3-2s)3(l + 2s)' 



(105c) 



These states above are exemplified in Figs. [T] and O For states such as P3;o,2: P4;0,2i and P4;o,3, convexifications (i.e. 
convex hull constructions) are needed; see Figs. [T] and [31 



D. Multi-qudits 

Now we consider multipartite qudit systems, in particular, the family of mixed states 

P^{Pk})^ E Pi\S{n-k)){S{n-k)\- 

The closest separable state that is in the symmetric subspace is of the form 

^* =E^fc|5(n;fc))(5(n;fc)|, 



(106) 



(107) 



and it can be constructed from the convex hull of product states 



(108) 
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FIG. 3: Comparison of F (solid), its convex hull, i.e., (partly dashed line, when convexification is necessary), and the 
numerical value of En for the states p4;o,2(s) (top) and p4;o,3(s) (bottom). 

Let us now describe how it is constructed and how the formula for the REE is obtained in a similar to the qubit case. 



Construct a separable mixed state by averaging over all phases, 



a{u)= rdcp,... r"#d|$(w,^))($(u,,?)|= VC5wt^...w^^|5(n;fc))(5(n;fc)|, (110) 
Jo Jo - 

k 

where, for the sake of convenience, we have defined (72 = n\/ki\ . . . kal. The relative entropy between p and a becomes 



Pi- 



k fe 1 « 



Minimizing this w.r.t. u's with the constraint J^i'^i = 1, we arrive at the solution 



(111) 



(112) 



This leads us to define the function 



PiPk) = J2Pk'^°S 



Pk 



(113) 



The relative entropy of entanglement for p(pg) is then conjectured to be the convex hull of the above expression. 
Theorem 2: 



MpiPk)) = ^^nPk)^ 



(114) 



and the closest separable state will be the corresponding convex hull of g{u). Whenever the function F is convex, 
there is no need for the last convexification procedure, and E'r = F. 
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Let us go on to prove this. As cr* in Eq. (|107p is constructed as the closest separable state to p{{pj:}) when restricted 
in the symmetric subspace, we have in particular that 



5/ 
dx 



(0,|<f.)(<i>.|)>0, 



for any \^s) of the form 



(115) 



k 



1=1 



This means that 



r *; 



(116) 



(117) 



for all qj > and J2i=i Ij — 1- Similar to the qubit case, if we can show that 



dl 

dx 



(o,l<f)(<f|)>o, 



for arbitrary product state |$), 



1=1 



for Ij.i = 1- We thus need to evaluate 



Note that 



niGS'„ 



fc2 



This gives 



(118) 



(119) 



(120) 



(121) 



(122) 



Taking the absolute value, we have 



/C2 



I < — X! (\/9n,(i),i ■ • ■ V9n,(fci),i) (V9n,(fci+i),2 • ■ • \/gnI(fcT+fe^) • • ■ {^J qii,{n-ka+i) .d ■ ■ ■ \/9n,(«),d) • (123) 



Recently, Carlen, Loss and Lieb [35| have shown an inequality regarding the permanent of a matrix. Using their 
results, the following inequality is easily seen to hold: 



3 E (V9ni(i),i • ■ ■ V9ni(fci),i) (V9ni(fci+i),2 ■ ■ ■ ^/Q^i{k^+k2),2) ■ ■ ■ {^/qni{n-ka+i),d ■ ■ ■ ^qni{n),d) < (qi) ■ ■ ■ (qd) , 
n, 

(124) 

where qi = J2j=i 'Zi,'/"-- With this, our conjecture is thus proved. 
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FIG. 4: Comparison of F (solid), its convex hull, i.e., Ek (dashed) for the state Pab{s) 



E. Examples 

Let us first at a three-qutrit example: n = 3, d = 4, and involving only two vectors a = (2, 0, 0, 1) and b = (1,1,1,0), 
to which the corresponding states are 



\a) EE i(|114> + |141) + |411)), 

\b) EE i(|123) + |132> + |213) + |231> + |312> + |321)). 



(125) 
(126) 



From Eq. ([4T|) their REE are 21og2(3) - 2 « 1.17 and 21og2(3) - 1 « 2.17, respectively. We shall now consider the 
mixture 



Pab{s) = s\a){d\ + il-s)\b)(b\. 
According to Eq. (|113p . the corresponding F function is 

9 

Fabis) = S l0g2 , ,2 + (1 - l0g2 
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(127) 



(128) 



(l + s)2 ' 2(1 -s2)' 

which is, however, not convex. This means that we need to construct its convex hull in order to obtain its REE, 

Er{s) = CO Fabis) = Fafc(l) + (1 - s){Fab{0) - Fab{l)). (129) 

This is shown in Fig. [D 

As another example, let us consider additionaly the state assocaited with c = (1, 0, 0, 2) defined via 



|c)^-(|441) + |414) + |144)). 

The state mixture of | a) and | c) , 

Pac{s) = s\a){a\ + {l-s)\c){c\, 
can be seen to possess REE equal to that of P3;i,2(s)- 



(130) 



(131) 



VI. .Biog < Br? 

Recall that for pure states we found the inequality Eiog < E^,- Does this inequality hold for mixed states? We do 
not know the complete answer to this question. From Inequality 6, we only have £'iog(p) — S{p) < £'r(p). However, 
we shall show that E'jog < -Er indeed holds for the two families of mixed states in Eqs. (p^l) and (|106p . as well as 
Eq. (fTiH . 
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A. Multi-qubits 



Let us begin with the multi-qubit states ((66)) . We first establish that the quantity £{{q}) that is used to obtain 
Slog is a lower bound on F{{q}), which is used to obtain Eyi- Recall that 0| 



£i{p}) = -2\og. 
By the concavity of log, we then have 



k n -^^n—k , 



cos y sm 



Hence 



21og2 



min —2 logn 



VP^ y^cos'^ B sin"-'= Q 



J2 Vp^ cos'' e sin"-'' e 



< 



^Pk l0g2 



Pk 



^ q'cos2'=6isin2("-'=)6l' 



- min^pfclogj 



Pk 



or equivalently 

f (M) < F(M). 

By using theorem 1 and by taking the convex hull of both sides of this inequality we have 

Eiog < En 

for the family of states ([55]) . 



(132a) 
(132b) 

(132c) 

(132d) 
(133) 



B. Multi-qudits 



Similarly, the relation i^iog < E'r also holds true for the multi-qudit states (|106p . and the proof is similar. 
Following the same idea in Ref . [2fl] , to calculate Ei^g for states (|106p , one first consider 



f({pJ) = -21og2 
By the concavity of log, we then have 



maxJ2VP^k\fcliV^iP ■ ■ ■ (v^)' 



-21og2 



k 



< 



Pk l0g2 



Pk 



Hence 



min —2 log2 



k k 



Pk 



<n^n^p,,log2— 



CSwi'^-i . . . Ud''" ' 

k 



or equivalently 

£{{Pf:}) < F{{p^}). 

By using theorem 1 and by taking the convex hull of both sides of this inequality we have 

Eiog < -Er 



for the family of states ([66]) . 

It would be interesting to know to what extent Eiog < -Er, holds. 



(134a) 



(134b) 



(134c) 

(134d) 
(135) 
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VII. SOME APPLICATIONS 



We know that the W state \W) = (|001) + |010) + |100))/\/3 is more robust than |GiJZ) ee (|000) + |111))/V2 agamst 
losing one of their constituent parties. This is also true for the particular family of n-qubit pure states {\S{n,k))}, 
the relative entropy of entanglement of which we have given in Eq. (|38p . Now, if we trace over one party we get a 
mixed {n — l)-qubit state: 



Tri|5(n,fc))(S'(n,fc)| 



n — k 



\S{n~l,k)){S{n~l,k)\ + -|S'(n-l,A:-l))(5(n-l,/c-l)|. 



(136) 



We have also given an expression for the relative entropy of entanglement for this mixed state. If we trace over m 
parties, the reduced mixed state would be a mixture of {|<S'(n — to, g))} [with q < (n — to)], and again we have its 
relative entropy of entanglement. For example, if we start with |S'(4, 1)), and trace over one party and then another, 
we get the sequence: 



|5(4,l))^P3;0,l(l/4)^P2;04(l/2), 



(137) 



for which we have given the corresponding relative entropies of entanglement in Eqs. ([55]) . (|104bp and (|103p . The 
entanglement at each stage is 



31og2 ( 2 



log2 
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9-31/4 



^log2 



or numerically. 



1.24511 ^ 0.433834 0.122556. 



(138) 



(139) 



Unlike n-GHZ states, these states can still remain entangled even if some of the qubits are lost along the way. 

Plcnio and Vedral 0] have derived a lower bound on the REE of a tripartite pure state pabc — in terms of 

the the entropies and REE of the reduced states of two parties: 



max{£;R(pAB) + S{pab), E-r{p Ac) + S'(pAc), -Br,(pbc) + S{pBc)} < Eji{pabc), 



(140) 



where pab = Trc(pABc) (and similarly for pAC and pec) and S{p) = — Trplog2 p is the von Neumann entropy. They 
have further found that this lower bound is saturated by |GHZ) and |W). Based on Theorem 1, we can actually show 
that for pi2...n = \S{n,k)){S{n,k)\ the inequality 



(141) 



is saturated, where p^^ i n = Tri(pi2...n) is the reduced density matrix obtained from pi2...n by tracing out the i-th 
party. The proof is as follows. As \S{n,k)) is permutation-invariant, there is no need to maximize over all parties, 
and we can simply take i = 1, obtaining the reduced state pn-i-k-i,k{k/n) as in Eq. (|136|) . As the corresponding 
function Fn-i-^k-i.k{s) of pn-i-^k-i,k{s) is convex for s G [0,1], we immediately obtain from Theorem 1 that, for 
Pn-i-k-i,k{k/n), 



E^Xpn-i-k-i,k{k/n)) = -log2 



l — k 



n—k 



k , k 
- l0g2 - 

n n 



n—k n—k 
log2 



= Sr. {\S{n,k))) - S {pn-i-k-i,k{k/n)) . 



(142a) 
(142b) 



Therefore, the bound in Eq. (|14ip is saturated for pi2...n = \S{n,k)){S{n,k)\. 

We remark that there are other applications that our results can be useful, such as (i) providing bounds on state 
discrimination by separable operations and (ii) constructing optimal entanglement witness, which have been discussed 
by Hayashi et al. [27| . and we refer readers to their paper. 



VIII. ANOTHER FAMILY OF MULTI-QUBIT MIXED STATES: DUR'S STATES 
Diir [l^l found that for TV > 4 the following state is bound entangled: 



Pn 



N 



(143) 
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where \^g) = (|0®^) + e^""|l'»^))/V2 is a iV-partite GHZ state; Pk = \uk){uk\ is a projector onto the state 
\uk} = |0)i|0)2...|l)fe...|0)jv; and = \vk){vk\ projects onto \vk) = |l)i|l)2 ■ • ■ |0)fc ■ • ■ |l)jv. For iV > 8 this state 
violates the (two-setting) Mermin-Klyshko-Bell inequahty [221 ; violation was pushed down to -/V > 7 by Kaszlikowski 
et al. [1^ for a three-setting Bell inequality; it was pushed further down to > 6 by Sen et al. ]40| for a functional 
Bell inequality. The phase aN in |5'g) can be eliminated by local unitary transformations, and hence we shall take 
a AT = in the following discussion. 

This state has been extended to the following family by Wei et al. [2lj |. 

PAr(x)EEx|vI/G)(vI/G| + ^^(P,-|-Pfe)' (144) 

fc=l 

and they found that for > 4 the state is bound entangled if < a; < 1/ {N +1) and is still entangled but not bound 
entangled if 1/{N -f 1) < a; < 1- They calculated the negativity and the geometric measure for this family of states, 
e-g-, 

2 

Eiog{pN{x)) = log2 j—^, (145) 

and conjectured their relative entropy of entanglement to be Eji{x) = x, for A'^ > 4, with one closest separable mixed 
state being 

a*{x) = 1{\^g){^g\ + I*g>(^gI) + ^ E (^'o + (146) 

k=l 

(147) 

We prove their conjecture here. As before we define the quantity 

f{z,as) = S{pn{x)\\{1 ~ z)a*{x) + zas), (148) 

where as is any separable state. In order to show that p*{x) is indeed the closest separable state to pm{x), it is 
sufficient to show that 

a/(0,^$)($|) ^ 1 _^°°^^Tr[(a* +t)-Vjv(^* +t)-»(a>|] >0, (149) 

where |$) is any separable pure state, which can be parametrized as follows, 

1$) = ®«^i(cos0,|O), -t-sin0,e'^^|l),), (150) 
and, without loss of generality, we can restrict ourselves to cos9j, sin^j > 0. By direct calculation, wc have 

r°° ^ _ 

/ dt((T*+<)-Vjv(f^*+t)"' = 2|*G)(*G|+V(Pfc+A), (151) 
Jo 



where 



fc=i 



and hence 

dz 

is valid for TV > 4 , as (see Ref. [IH for the proof) 



1-.9A'>0 (152) 



gN = (ci • • • Cat + ,Si • • • Sat 



2 



N 



'^{{ci---Sk - ■■cn)'^ + {si---Ck--- snY] < 1, (153) 



k=l 



where we have simplified the notation by using Ci = cos 6i and Si = sin 6i. Therefore, (7*{x) is indeed the closest 
separable state to pn{x) and hence Er(^pm{x)) = a; for > 4. 

We remark that for this family of states, the relation E\og < E^ holds, as log2 (2/ (2 — x)) < x for < a; < 1. 
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IX. CONCLUDING REMARKS 



We have proved conjectures on the relative entropy of entanglement (REE) for two families of multipartite qubit 
states. Thus, analytic expressions of REE for these families of states can be straightforwardly obtained. The first 
family of states we consider are composed of mixture of some permutation-invariant multi-qubit states. We have 
also generalized the results generalized to permutation-invariant multi-qudit states, and have given the expression for 
their relative entropy of entanglement. The second family of states contain Diir's multipartite bound entangled states. 
Along the way, we have reviewed inequalities connecting the relative entropy of entanglement to the robustness of 
entanglement and the geometric measure of entanglement, and have slightly extended previous discussions. 

It is possible that our results on the relative entropy of entanglement can applied to the checking of the consistency 
of some equalities and inequalities 0, regarding minimal reversible entanglement generating sets (MREGSs). 

These equalities and inequalities concerning MREGS usually involve only the von Neumann entropy and the regu- 
larized (i.e. asymptotic) relative entropy of entanglement of the pure state and its reduced density matrices. The 
results we have in the present Paper concern only the non-regularized version of the relative entropy of entanglement, 
and hence, can only reach weaker conclusion. Therefore, a major challenge is to extend the ideas contained in the 
present Paper to the considerations of the regularized version of the relative entropy of entanglement. Moreover, 
knowledge of the regularized version is also important for providing bounds on the random bipartite entanglement 
recently investigated by Fortescue and Lo [43| . 

In Ref. [iBl the geometric measure of entanglement was known for not only the family of mixed states 
J2kPk\S{n,k)){S{n,k)\ but also the family of pure states (formed by superposition) J2k^k\S{n,k)). For the rel- 
ative entropy of entanglement, we still do not yet have analytic expressions for the latter, except for the basis states 
\S{n,k)). It will be desirable to consider how we can obtain the relative entropy of entanglement for this family of 
pure multi-qubit states, as well as the corresponding family of pure multi-qudit states '^i:ctf:\S{n,k)). 

Indeed, investigating the relative entropy of entanglement for states in this family '^j:cti:\S(n,k)) can serve as a 
first step towards obtaining the regularized version of the relative entropy of entanglement for states such as \S{n, k)) 
and \S{n, k)), as we now illustrate. Let us consider two copies of the W state (shared among parties A, B, and C), 

\W^^) = \W)a,b,c, ® \W)a,b,c,, (154) 

where 

|W) = i=(|001) + 1010) + 1100)). (155) 
We can expand and re-label the states as follows 



|0)ai®|0)a. ^ \1)a, (156a) 

\0)a,®\1)a, \2)a, (156b) 

\1)a,®\0)a, ^ |3)a, (156c) 

\1)a,®\1)a, ^ |4)a, (156d) 

and similarly for parties B and C. This gives us 

IW^') = ^\a)ABC + ^\b)ABC, (157) 

where |a) and |6) are defined in Eqs. (|125p and p26p . respectively. Thus, the two copies of the W-state is identical 
to a superposition of two four-level three-party states, exactly of the form ag|5'(n, k)). In general, a state with m 
copies of, say, \S{n, k)) is also of this form, and hence, the knowledge of the relative entropy of entanglement for the 
former enables the knowledge of the regularized relative entropy of entanglement for the latter. 

In discussing the symmetry of the separable states, we have essentially characterized the family of states invariant 
under such symmetry (Pi and P2). It would be desirable to quantify the entanglement of such family. This family 
includes states that are generalizations of Werner states, which invariant under average action of U (d)®" in multipartite 
systems. An interesting question arises: can we characterize the distillability of multipartite states just as in bipartite 
case, where the question can be reduced to that for Werner states? 
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